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1. Motivation 

The phenomenology of QCD is governed by two prominent features, confinement and sponta- 
neous breaking of ciiiral symmetry. As one increases the temperature T above some critical value 
Tc, the theory becomes deconfined and chiral symmetry is restored. This suggests that there could 
be a relation between the two phenomena. Establishing or ruling out such a relation would be a 
major insight into key mechanisms of QCD. 

The finite temperature transition of pure gauge theory, where the system changes from the 
confined {T < Tc) into the deconfined phase {T > Tc), can be understood as spontaneous breaking 
of the center symmetry [[I]], and the Polyakov loop is a suitable order parameter with (P) = 
for T < Tc and (P) 7^ for P > r^. (see Fig. [I]). Concerning chiral symmetry breaking the order 
parameter is given by the chiral condensate with (yy) / for P < P^ and (v^V^) = for P > P^. 
The chiral condensate in turn is related to the spectral density of the Dirac operator by the Banks- 
Casher formula [^]: (yy) = ~^P(0)' where p(0) is the spectral density at the origin. When 
increasing the temperature above the critical value, the spectral density at the origin vanishes and 
with it the chiral condensate (see Fig. ^). 

In a series of recent papers, |^-[^, the Polyakov loop, or other quantities that serve as order 
parameters for the breaking of the center symmetry, were related to spectral sums of differen- 
tial operators on the lattice, in particular the Dirac and the covariant Laplace operators. In this 
way also confinement is related to spectral quantities, as is chiral symmetry breaking through the 
Banks Casher formula. Since eigenvalues can be divided into IR and UV regions in a natural way, 
the spectral representations allow one to analyze whether confinement is dominated by IR or UV 
modes. In this paper we review and extend our contributions to this enterprise. 

2. Derivation of spectral sums 
2.1 Preliminaries 

The thin Polyakov loop averaged over space is given by 



where L(x) is a straight line of links in temporal direction closed around the periodic boundary and 
Trc denotes the color trace. Obviously, this quantity has a rather singular support. Due to this fact 
it is known to have a poor continuum limit with large renormalization effects. 

To address this problem, below we will define a new observable P^'^^ which we refer to as the 
dressed Polyakov loop with winding number q. It is a collection of loops all with the same number q 
of windings around compactified time. So, for example, P^'' is a sum of the thin Polyakov loop and 
all other loops winding once around the lattice in time direction, where the latter loops may have 
arbitrary complicated detours in spatial directions and are weighted proportional to their length. 

In the following two sections we derive the spectral sums for thin and dressed Polyakov loops. 
For this purpose it is convenient to summarize our conventions: The staggered Dirac operator is 
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Figure 1: Scatter plot for the thin Polyakov loop in the complex plane for an ensemble of quenched SU(3) 
configurations. As the temperature increases (left to right) the Polyakov loop acquires a non-vanishing 
expectation value. 
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Figure 2: Histograms of the spectral density of the staggered lattice Dirac operator as a function of the 
modulus of the eigenvalues. Below the critical temperature the spectral density at the origin is non-zero. 
Near Tc the density at the origin decreases and vanishes above the critical temperature. 



given by 

1 

D{x,y) = md,^y + X £ % W [U^ W4+A,y -U^ix-jX) ^4--A,y] ' (2.2) 

where 77^(x) = (—lyfi+'^z+.-.+JCu-i jjjg staggered phase. The covariant Laplace operator on the 
lattice is defined as 

4 

A{x,y) = (8 + m2)5,,y- £ [[/^W5,+^,, + [/^(x-/i)^5,_^,,] • (2.3) 

Ai=l 

Below we will use the staggered Dirac operator for representing the thin Polyakov loop, while the 
dressed loops are constructed from the Laplace operator. We stress, however, that it is an easy 
exercise to generalize our formulas, such that both the thin and the dressed Polyakov loops can be 
obtained from both differential operators. The derivation can be extended further to arbitrary lattice 
operators with only nearest neighbor interaction (such as Domain Wall Fermions). The advantage 
of these operators is that they are numerically cheap. 
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2.2 Thin Polyakov loops 

In order to connect the thin Polyakov loop to spectral sums we follow the derivation in [^]. We 
consider powers of the staggered Dirac operator (^^. The n-th power of the operator then contains 
paths up to length n, dressed with a product of n links. For the special case where we consider the 
operator to the power Nt, where Nt is the time extent of our lattice, we obtain 
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N, Y N,-l 

Tr,D'''{x,x) = -±^Tr,YlU4{x,s)--^Tr,llU4{x,Nr-s) 

^ ' s=\ ^ ' .v=0 

+ Other loops, trivially closed 
2^ [^(■?) — ^*(-^) ] + other loops, trivially closed. 



(2.4) 



(2.5) 



From the set of all loops that contribute we have singled out the thin Polyakov loop and its complex 
conjugate. These two straight loops are the only ones that can close around compactified time. All 
others are 'trivially closed', i.e., they do not wrap around the lattice in temporal direction. 
The key observation is that a change of the temporal boundary conditions 



U^ix.Nt) ^ zt/4(x,A^,) , \z 
where z is a complex phase, only affects the Polyakov loops 
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(2.6) 



(2.7) 



while the trivially closed loops remain unchanged. After such a change the expression ( |2.5[ ) reads 



Tr,lf^'{x,x) = ■:^[zL{x) -z*L*{x)] + other loops, 



(2.8) 



where is the Dirac operator on a configuration where the temporal boundary conditions are 
changed by a factor z- Averaging this expression over space and time we obtain 
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Yl.{zP-z*P*+X), 



(2.9) 



where X is the sum of all trivially closed paths. We make use of this behavior to cancel the un- 
wanted trivial contributions and project onto the thin Polyakov loop by taking linear combinations 
(coefficients at) of this expression for three different boundary conditions zi,Z2 and 23 
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One can show that this equation has a unique solution as long as the boundary conditions Z[,Z2 and 
Z3 are different. For the particular choice zi = 1 and 22=^3= = z the thin Polyakov loop is 
then given by 



2N, 

P = — 
V4 



N, 



(0 



N, 



(2.12) 



where the Xy'^ are the eigenvalues of the staggered Dirac operator for a given boundary condition 
z. Similar expressions can be obtained for other differential operators on the lattice. 

2.3 Dressed Polyakov loops 

In order to obtain the spectral sum for the dressed Polyakov loops we now use the Laplace operator 
and apply a different strategy than used for the thin loops. First, we rewrite the covariant Laplace 
operator on the lattice from Eq. (2.3) in the following way 



4 , 



-[1-kH] , H{x,y) 



I 

il=\ 



(2.13) 



where the hopping matrix H collects all terms which connect nearest neighbors on the lattice and 
the hopping parameter K is related to the bare mass m via K = (8 + m^)^^ The inverse Laplace 
operator A^' {x,x) can be expressed as geometric series in terms of powers of H, 



A ^(x,x) 



(2.14) 



Since the hopping matrix H contains only terms that connect nearest neighbors, the powers H^{x,x) 
in ( 2.14| ) correspond to closed loops of length j, again dressed with the corresponding products of 
the link variables. Thus we can organize A^' in terms of loops and order these with respect to their 
winding around time. We obtain (after taking the trace) 



TrA 



-1 



£Tr,A-Hx,x) = £ £ ^I'l+'Tr, nf/.(3') 



(2.15) 



where is the set of loops that wind exactly n-times around the time direction and |/| denotes 
the length of a loop /. It is important that now the sum on the right hand side contains loops 
of arbitrary length, with all possible numbers of windings. To disentangle the different winding 
numbers we introduce U(l)- valued temporal boundary conditions for the Laplace operator which 
are most conveniently implemented by the replacement 



U4{x,Nt) , 



(2.16) 



which corresponds to a parameterization z = e'f in ( P^ . Each loop acquires a phase factor e"f", 
corresponding to its winding number n: 



Tr Ac 



-1 



(2.17) 
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Integrating over cp with a factor of e projects to winding number q and we end up with the 
dressed Polyakov loop P^^) (which we normalize with V = N^Nt), 



-E 



1+1 



V 



(2.18) 



After the projection the loops still can be arbitrarily long but all have the same winding number q. 
The individual loops are weighted with a factor K"l'l+' , where |/| is the length of the loop. This leads 
to an exponential suppression of long loops and the rate of the suppression is determined by the 



value of K. Nevertheless arbitrary long loops may contribute and the formula ( |2. 18| ) is a sensible 
definition of dressed loops. 

The final step to make the connection to spectral properties of the staggered Dirac operator is 
to write the trace over A^^ as a spectral sum giving 
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(2.19) 



where A j'''^ is the j'-th eigenvalue of the hopping matrix H ( |2. 13| ) for boundary condition angle <p 



The dressed Polyakov loop transforms in the same way as the thin Polyakov loop. Thus the dressed 
Polyakov loops are proper order parameters for the breaking of the center symmetry. They are 
also not as singular as the thin Polyakov loops and thus are expected to have a better continuum 
limit with smaller renormalization effects. In an analogous way dressed Polyakov loops can be 
constructed also with the eigenvalues of the Dirac operator instead of the Laplace spectrum. 

We stress that both equations ( 2.12 ) and ( 2. 19| ) are exact results that hold for individual gauge 
configurations. However, for a practical evaluation with numerically generated gauge configura- 



tions the two formulas have a different status. Equation ( |2.12| ) for the thin loops holds to ma- 
chine precision when the result of the spectral sum is compared to the direct evaluation of the thin 



Polyakov loop from the gluonic definition ( |2.1[ ). The expression ( |2.19[ ) for the dressed loop on the 
other hand contains an integral over a continuum of boundary conditions. That this integral can 
indeed be approximated numerically in a sensible way will be demonstrated in the next section. 
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12^ X 6 


0.098(1) 
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Table 1: Ensembles used in our simulations. For the 6^ x 4 lattices 2000 configurations are available, for 
12^ X 4 and 12^ x 6 we have used 20 configurations. 
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3. Numerical analysis 
3.1 Simulation details 



For the numerical analysis of our spectral sums ( 2.12 ) and ( [2.19 ) we use quenched SU{3) gauge 



configurations generated with the Liischer-Weisz action 0. We work on lattices with different 
sizes ranging from 6^ x 4 up to 12^ x 6 at temperatures below and above the critical value (see 
Table |l]). For those we compute complete spectra of both the staggered Dirac operator and the 
hopping matrix H of the Laplace operator for up to 16 different boundary condition. The lattice 
spacing was set [§] with the Sommer parameter ro = 0.5 fm. All error bars we show are statistical 
errors determined with single elimination jackknife. 

3.2 Thin Polyakov loops 

For the thin Polyakov loops we restrict ourselves to the three boundary conditions zi = 1 and 
Z2 = Z3 = e'^ = z which ai^e used to derive the spectral sum ( 2.12| ). A different choice would 
change our results only marginally. 

After computing all eigenvalues of the staggered Dirac operators we order them w.r.t. to their 
absolute value and organize them in bins. For each bin we calculate a number of observables. 
The first observable is the distribution of the eigenvalues, which we show in Fig. |3| for different 
spatial and temporal extents of the lattice and for three different temperatures above, below and 
approximately at the critical value T^. The l.h.s. set of those plots shows the distribution as a 
function of the size of the eigenvalues in lattice units, while the r.h.s. plots are in MeV. One can 
clearly see that the spectral density at the origin vanishes and chiral symmetry is restored as one 
increases the temperature above the critical value. The density of eigenvalues reaches a maximum 
after about three quarters of the eigenvalues and then quickly drops towards the UV cutoff. 

Fig. ^ is particularly important when one analyzes which part of the spectrum contributes to 



the spectral sum ( |2. 12D for the thin Polyakov loop. For the interpretation of such an analysis it is 
necessary to take into account the global distribution of the eigenvalues, and this is exactly what 
Fig. m provides. 



In the spectral sum ( [2.12| ) the thin Polyakov loop emerges through a shift of the eigenvalues as 
the boundary condition is changed. It is interesting to analyze how strongly different eigenvalues 
of the spectrum are shifted when the boundary condition is changed. We quantify this question by 
studying the averaged shift of the eigenvalues given by 

s{X) = + + |A,-A,.|)/3 . (3.1) 

Since above the phase transition the thin Polyakov loop obtains a finite value while in the confined 
phase it is approximately zero, the shift of the eigenvalues should change as one increases the 
temperature. This is exactly what we find in our results for the average shift s{X) presented in 
Fig. The individual data points are obtained by averaging over all eigenvalues in a bin and then 
over all configurations in the ensemble. 

It is obvious, that over the whole range of the spectrum the average shift is almost zero below 
the critical temperature. For high temperatures, however, we find large shifts for low-lying eigen- 
values and little or no shift when the UV region is approached. The small shift at the UV end of 
the spectrum disappears when finer and larger lattices are used. 
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Figure 3: Distribution of the eigenvalues of the staggered Dirac operator for different lattices and different 
temperatures. The l.h.s. plots are in lattice units, while on the r.h.s. we use MeV. 



After understanding the global distribution of the eigenvalues and their shift under a change 
of the boundary conditions, we can start to analyze the individual contributions to the spectral sum 



(2.12). This individual contribution is given by 



2N, 



(3.2) 



Our results are shown in Fig. |5|, where we plot the absolute value of the individual contri- 
butions after normalizing them with P, i.e., we divide by the value of the thin Polyakov loop. It 
is obvious that mainly the UV modes contribute to the thin Polyakov loop. This observation is 
essentially independent of the temperature and the size of the lattice. This is a finding which is not 
a-priori obvious, given the fact that it is the IR modes that show the largest shifts in Fig. ^. On the 



other hand, the large power Nt in the spectral sum ( |2.12| ) drastically enhances the contributions of 
the UV eigenvalues. 

Interesting are also the dips which form with increasing temperature at about 1.5a for Nt =4 
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Figure 4: Average shift of the eigenvalues calculated as defined in Eq. (3.1) as a function of X. Again we 
use lattice units in the l.h.s. set of plots and MeV on the r.h.s. 



and at 0.8a and 1.9a for Nt = 6. It can be shown that these dips are correlated with a change of the 
direction of the shift of the eigenvalues for z- valued boundary conditions relative to the eigenvalues 
computed with periodic boundary conditions. 

The plot in Fig. || shows that the individual contributions are biggest at the UV end of the 
spectrum. On the other hand, Fig. |3| shows that for the largest eigenvalues the density shows a 
sharp drop. It is an interesting question which effect, size of contribution versus density, wins out. 
This question can be addressed by considering only partial sums of ^A% . In Fig. H we show 
such partial sums as a function of the cutoff, i.e., as a function of the largest eigenvalue included in 
( 2.12 ). Again we divide by the Polyakov loop P, such that when all eigenvalues are summed, i.e., 
at the largest value of A in Fig. ^, the curve approaches 1. 

For the 12^ x 6 lattice, we find that for low and medium temperatures the partial sums over- 
shoot by nearly a factor of two before they reach the correct value at the largest values of A . Thus 
the behavior of the partial sums is non-monotonic. On the other hand, for small temperatures the 
Polyakov loop should vanish. Thus the overshooting is not problematic, since the accumulated 
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0,0 0,5 



Figure 5: Contribution (3.2) of individual eigenvalues to the spectral sum for the thin Polyakov loop given 
by Eq. (2.12). As before we use lattice units in the l.h.s. set of plots and MeV on the rh.s. 



contributions stay approximately zero over the whole range of the spectrum. We also find that the 
thin Polyakov loop is recovered only after all UV modes are included in the spectral sum. Again 
we see the appearance of dips at exactly the same values which we have already identified in Fig. ^. 

So far we have only considered absolute values of quantities. However, the Polyakov loop and 
also the spectral sums are complex numbers and we can analyze their relative phase. In Fig. we 
plot the angle between the truncated sums and the resulting thin Polyakov loop as a function 
of the cutoff used in the truncated sum. We show data for a 12^ x 4 and a 12^ x 6 lattice, both 
at r > Tc. We find that on the former lattice the accumulated contributions start with the wrong 
sign and recover the correct sign at the end. On the other lattice, although starting and ending with 
the correct sign, there is an intermediate region in which the truncated sum has the wrong sign. 
A comparison with lattices of different spatial size but the same temporal extent reveals, that this 
behavior depends only on the temporal extent. 

A possible explanation for this phenomenon is the fact that the eigenvalues A of the staggered 
Dirac operator are purely imaginary. Thus A^" is a real positive number, while X^"^^ is real but 
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Figure 6: Partial sums of the spectral formula (2. 12), plotted as a function of the cutoff. 



negative. If we consider only the first term in the truncated sum we find: 



'yN, 



N, 



N, 



N,- 



2N, 

— [a + {z*+z)&] 



2^' _ a-t~a-=a 

— [a + z a^+za,.] ~ =- 

2^' ^(<0N,=4n asa,a>0 

— a — a < , 
V4 ^\>0Nt=4n + 2 asa,a<0 



(3.3) 



where we assume that the boundary conditions with phases z and z* lead to approximately the same 
lowest eigenvalue a, which is larger in size than the eigenvalue at periodic boundary conditions 
^. Interestingly, this last relation between the relative sizes of the eigenvalues for non-trivial and 
periodic boundary conditions seems to change in certain regions in the spectrum of the staggered 
Dirac operator, thus creating the observed sign change in the truncated sums. The points where 
these sign changes occur depend only on the time extend of the lattice. 
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Figure 7: Phase shift of the truncated sums relative to the thin Polyakov loop for T > T^-. The upper plots 
are for a 12^ x 4 lattice, while the lower plots are for 12^ x 6. 



3.3 Dressed Polyakov loops 

For the thin Polyakov loops the numerical analysis of the last section demonstrates that the 
spectral sums for the thin loops are predominantly built up from the UV modes. However, as we 
have already speculated above, part of the UV dominance might come from the singular nature 
of the support of the thin loop. Thus it is interesting to perform the same spectral analysis for 
the dressed Polyakov loop. Such an analysis is planned for a future publication. Here we focus 



on the feasibility of computing the integral over the boundary angle in ( [2. 19[ ) numerically and 
establish that the dressed Polyakov loop is indeed an order parameter for the breaking of the center 
symmetry. 

For the dressed Polyakov loops we compute complete spectra of H for 16 values of the bound- 
ary angle (p £ [0, tt). The integrand in the spectral sum is given by (omitting the normaliza- 
tion and the phase exp{—i(pq)) 



Si(p) 



1 3^ 
^ k 1 



(9) ■ 



(3.4) 



Integrating S{(p) over <p with the phase exp(— /<p^) projects to the loops with winding number q. 
The important technical question is how smooth a function of (p this integrand is. 

Fig. ^ shows the integrand S{(p) averaged over 20 configurations above Tc with real Polyakov 
loops (l.h.s. plots) and 20 configurations with complex Polyakov loops (r.h.s.). The individual 
curves are for different values of the mass parameter m in the vicinity of the critical mass where 
the Laplace operator ceases to be invertible. In the bottom plots a constant was subtracted. 
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Figure 8: The spectral sum S{(p) for the Laplace operator as a function of (p. The l.h.s. plots are from 
an average over 20 configurations with real Polyakov loop, the r.h.s. plots are for 20 configurations with a 
Polyakov loop with phase ^ exp(/2;r/3). In the top plots the full sum 5'(^) is shown, while in the bottom 
plots the constant c = / d(pSL{(p) was subtracted. 



The plots show a clear cosine behavior for the ensemble with real Polyakov loops (l.h.s.) and a 
shifted cosine for the complex ensemble (r.h.s.). For low-lying modes a similar dependence on the 
phase at the boundary was observed in |[T^]. It is obvious that dividing the interval [0,7r) into 16 
sub-intervals should be sufficient for a reliable numerical estimate of the (p-integral in the case of 
q = I, where the phase factor is simply coscp — / sin (p. This is the most interesting case of dressed 
Polyakov loops with a single winding. 

Fig. |]estabhshes that a numerical integration over <p is feasible. We implement this integration 
for the case of ^ = 1 using Simpson's rule. The resulting dressed Polyakov loop P^^^ transforms 
under center transformations in the same way as the thin Polyakov loop P defined in (p|). Thus it 
should be an order parameter for breaking of center symmetry, but one even expects that it shows 
a similar behavior as shown in Fig. |l| for the thin loop. In order to test this expectation we use 20 
configurations below Tc and 20 above Tc. The latter all have a Polyakov loop which is essentially 
real, i.e., configurations in the right leg of the distribution shown in the r.h.s. plot of Fig. |l]. On 
these two sets of configurations we compare the thin loop with the dressed loop at two different 
values of k. The corresponding results in the complex plane are shown in Fig. |^ 

The plot shows clearly that for T < Tc both, the thin and the dressed Polyakov loops cluster 
around the origin, while above T^ they develop a non-vanishing value near the positive real axis. 
Even the pattern in the distribution of the data points is similar for thin and dressed loops. 
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Figure 9: Scatter plots of thin (l.h.s. plot) and unsealed dressed Polyakov loops (center and l.h.s. plots) 
below (blue symbols) and above Tc (red symbols) for 20 configurations on a 12^ x 6 lattice. 



The analysis of Fig. ^ shows that the dressed Polyakov loop ( 2.19| ) is a suitable order parameter 
for the breaking of center symmetry. It is less singular and is expected to have a proper continuum 
limit. We are currently analyzing the volume and temperature dependence of the dressed loops and 
will then return to studying which part of the spectrum predominantly contributes to the spectral 
sum for the dressed loops ( |2.19 ). 
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